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Abstract 

Newton’s third law states that any action is countered by a reaction 
of equal magnitude but opposite direction. The total force in a system 
not affected by external forces is thus zero. However, according to the 
principles of relativity a signal can not propagate at speeds exceeding 
the speed of light. Hence the action cannot be generated at the same 
time with the reaction due to the relativity of simultaneity, thus the 
total force cannot be null at a given time. The following is a contin¬ 
uation of a previous paper [T] in which we analyzed the relativistic 
effects in a system of two current conducting loops. Here the analysis 
is repeated but one of the loops is replaced by a permanent magnet. 
It should be emphasized that although momentum can be created in 
the material part of the system as described in the following work 
momentum can not be created in the physical system, hence for any 
momentum that is acquired by matter an opposite momentum is at¬ 
tributed to the electromagnetic field. 
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1 Introduction 

Among the major achievements of Sir Isaac Newton is the formulation of 
Newton’s third law stating that any action is countered by a reaction of 
equal magnitude but opposite direction [a [3]. The total force in a system not 
affected by external forces is thus zero. This law has numerous experimental 
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verifications and seems to be one of the corner stones of physics. However, 
by the middle of the nineteenth century Maxwell has formulated the laws 
of electromagnetism in his famous four partial differential equations I31E11I1 
which were formulated in their current form by Oliver Heaviside [8]. One of 
the consequences of these equations is that an electromagnetic signal cannot 
travel at speeds exceeding that of light. This was later used by Albert 
Einstein [H [71 E] (among other things) to formulate his special theory of 
relativity which postulates that the speed of light is the maximal allowed 
velocity in nature. According to the principles of relativity no signal (even 
if not electromagnetic) can propagate at superluminal velocities. Hence an 
action and its reaction cannot be generated at the same time because of the 
relativity of simultaneity. Thus the total force cannot be null at a given 
time. In consequence, by not holding rigorously the simultaneity of action 
and reaction Newton’s third law cannot hold in exact form but only as an 
approximation. Moreover, the total force within a system that is not acted 
upon by an external force would not be rigorously null since the action and 
reactions are not able to balance each other and the total force on a system 
which is not affected by an external force in not null in an exact sense. 

Most locomotive systems of today are based on open systems. A rocket 
sheds exhaust gas to propel itself, a speeding bullet generates recoil. A car 
pushes the road with the same force that is used to accelerate it, the same 
is true regarding the interaction of a plane with air and of a ship with wa¬ 
ter. However, the above relativistic considerations suggest’s a new type of 
motor in which the open system is not composed of two material bodies 
but of a material body and held. Ignoring the held a naive observer will 
see the material body gaining momentum created out of nothing, however, 
a knowledgeable observer will understand that the opposite amount of mo¬ 
mentum is obtained by the held. Indeed Noether’s theorem dictates that 
any system possessing translational symmetry will conserve momentum and 
the total physical system containing matter and held is indeed symmetrical 
under translations, while every sub-system (either matter or held) is not. 

In this paper we will use Jehmenkos equation [6] discuss the force be¬ 
tween a current loop and a permanent magnet. In this respect the current 
paper differs from a previous one [T] which discussed the case of two cur¬ 
rent carrying coils. The current conhguration may seem attractive since a 
permanent magnet does not require a power source. 
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2 The Magnetic Current Density 

Let us begin by writing down Ampere’s law for the magnetic field H : 

V X H = J + dtD. (1) 

We use the standard notations of vector analysis, J is the electric current 
and D is the displacement held. In matter the magnetic held is related to 
the magnetic hux density B as: 


H = — -M. (2) 

/^o 

In which M is the magnetization and fio is the vacuum permeability. Hence; 

V X {—)-V X M = J + dtD. (3) 

fJ'O 

We now dehne the magnetic current as: 

xM. (4) 

To obtain: 

V X H + </ + dtD). (5) 

Thus, the magnetic current density Jm plays the same rule as the current 
density J. 


3 Configuration 

We will assume a uniform magnetization M = Mz{r,z)z conhned between 
the planes z = —hi and z = —h such that h < hi and hi — h = L. The 
magnet has a cylindrical shape such that the magnetization is conhned to 
as radius r < i? 2 - So that: 


Mz = 


Mo — hi < z < —h 
0 otherwise. 


&: 0 < r < i? 2 - 


In terms of the step function u{x) 


( 6 ) 


u{x) = 


1 X >0. 

0 X < 0. 
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(7) 





Figure 1: A cylindrical magnet (blue) and a current loop (red) above it. 
(three different sections) 


the magnetization can be written as: 

M 2 = Mo {u{z + hi) - u{z + h)) u(i ?2 - r). (8) 

The magnet and a current loop above it are depicted in figure [H 

We can now calculate the magnetization current by evaluating equation 
dH using cylindrical coordinates: 

Jm = -[deiM^) - dz{rMo)]r+ [dz{Mr) - dr{M^)]6 
r 

+ ^[drirMe) - de{Mr)]z = ^d0{Mz)r - dr{Mz)9. (9) 

Since = Mz{r, z) we obtain: 

Jm = —drMzO. ( 10 ) 

Moreover, since the Dirac delta function 5{x) is related to the step function 
by the formulae: 

/ X 

5{x')dx'. (11) 

-CO 
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We obtain: 


^ ^ du{R2 - r) d{R2 - r) 

- ■■> = d{R, - r) dr = 

Using equation (fT^ and equation ([8]) in equation (fTOjl we obtain: 

Jm = Mq{u{z + hi) — u{z + h))6{R2 - r)9 
The above can also be written in the following form: 




rz+h 


Jm = Mo{ 5 {z')dz'- 6{z')dz')6{R2 - r)e 


—oo 

z-\-hl 


= Mo 


rz-\-n 

J z+h 


5{z')dz'5{R2 - r)e. 


Now consider a change of variables: 


z" = z' — z. 


This will lead to: 


JM{r,z) = Mo 


rhl 

/ 

Jh 


dz" = dz' 


z” + z)dz'' 5 {r - R2)6. 


Making an additional variable change: 


z'” = —z'' 


Results in: 


JMir,z) = Mo 


/-%( 


dz'” = -dz'‘ 


z - z'")dz'" 5 {r - R2)e. 


( 12 ) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


Hence we see that the magnetization current density is equivalent to a con¬ 
tinuum of loop currents each with a current density: 


f'^{r,z,z') = Mo6{z - z')6{r - R2)9. 


(19) 


We notice that the magnetization Mq replaces the current I appearing in 
the analog expression for charge current density. Hence we can write: 


^ r— j* 

JM{r,z)= J^^{r,z,z')dz'. 
J-hl 


( 20 ) 
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4 Force Calculations 


Since in a previous paper [T] we have calculated the relativistic total force 
acting on a system of two coils (see equation (44) of [T]) we will take ad¬ 
vantage of this calculation taking into account that the magnetic current 
density is static. For a single loop we have: 


^(1) ^ M 
‘ Svr c 


i2,2hh - -(-, 




( 21 ) 


In which 2 ) dehned in equation (39) of [T] and we have replaced I 2 
with Mq. The total relativistic force is an integration over the contribution 
of a continuum of current loops: 


Ft = 


^ — h 


'-hi 


OTT C 


r‘ — h 


-hi 


KS^dz' = ^{^lfMoiiKi2,2. ( 22 ) 

’ OTT C 


In which iFi 2,2 is an integration the K factor of all magnetization current 
loops. Now we wish to calculate 2 ) equation (39) of [T] 

which is: 

^S,2= Rdh ■ dh . (23) 

Notice that R = Xi — X 2 is a difference vector between the location vectors 
of current elements on the charge current and magnetization current loops. 
Hence using suitable variables: 

= (i?i cos 01, sin^i, 0), X 2 = (i ?2 cos 02) .R 2 sin02,2^)- (24) 

We obtain: 

R = Xi — X 2 = {Ri cos 01 — R 2 cos 02, i?i sin 0i — R 2 sin 02 , —z')- (25) 


And thus R is: 

R = \/{Ri cos 01 — R 2 cos 02)^ + {Ri sin 0i — R 2 sin 02)^ -|- z'‘^ 

= ^Rj + Rl- 2 R 1 R 2 cos (01 - 02 ) + z'2. (26) 

Now we can calculate the vector line elements: 


dll = i?id0i0"i = i?i(i0i(—sin 01 , cos 01 ,0), 

dh = i?2d0202 = .R2d02(-sin02,cos02,O). (27) 
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The scalar product of those line elements is: 

dll ■ dh = RiR2d6id92 cos( 0 i — 6 * 2 ). (28) 

Combining the above results we can write the K integral of equation ()23l) 
as: 

= r de 2 d9icos{9i-e2)R{ei,92,z'). (29) 

’ h Jo Jo 

We now make a change in variables: 

e' = 9^-92, d9' = d9i. (30) 


The second above equation is correct since 02 is constant for the 9i integral. 
Hence: 

K^^ 2 iz') = T" d92 r d9' cos {9')R{9' + 02 , 02 , z'). (31) 

’ d Jo J-g^ 

In the above: 


^ R (i?i cos (0'+ 02 ) — i ?2 cos 02 , i?i sin (0'+ 02 ) — i ?2 sin 02 , —z') 

~ R~ y/R:i +Ri- 2 R 1 R 2 cos 0' + z'2 

(32) 

Now we notice that all the function contained in the integrand of equation 
(|3T]l are periodic in 0 ' with a period of 27r hence we can replace d9' —>• 

d9'. The following step would to change the order of integrals performing 
the 02 integral hrst and noticing that all the functions which are periodic in 
02 have a null contribution to the integral. Hence we obtain: 



2ttRiR 2 d9'cos 9'z' 

Jo i/Rlf + Ri- 2 R 1 R 2 cos 0 ' + z '2 


(33) 


Now we sum up contributions to the K factor from all loops: 


i?i 2 , 2 = [ dz'K[l\iz') = r^d 0 'cos 0 ' [ ^ dz' 

j-hi ’ d Jo J-hi 


Va^~+ z‘ 


/2 


(34) 

in which we define: o? = + i?! ~ 2 R 1 R 2 cos 9'. A simple integration leads 

to: ^ 

Ki 2,2 = —^ J d9' cos 9'{-\/a'^ + — ija'^ + d1)z. (35) 
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Now we can write: 


\/+ R\ + R\ — 2 R 1 R 2 cos 9' 

= y/R^R^J——— +^ + ^ - 2cos 9'. (36) 

y K 1 K 2 1 x 2 -fti 

And define; 

/r 2 i?i i ?2 

“ R 1 R 2 R 2 Ri 

Hence: 

\/q ;2 + /l^ = i?ii? 2 V^ — 2 cos 0 ', 

In terms of the above definitions 

^ 2'k{RiR2Y-^ ^ 

-fi^i 2,2 =-- {g{b) - g{bi))z. (39) 

In which: 

r27r 

g{b) = / dO' cos OWb — 2 cos 6'. (40) 

Jo 

It is obvious that when 6 S> 2 the respective part of the integral vanishes 
(the same is true for bi) this is evident since the integral is performed over a 
constant time a cosine function of period 27r. It also clear that 6 is a sum of a 
factor dependent on the magnet vertical dimensions and a factor dependent 
the ratio between Ri and R 2 which we denote s = ^ . The second factor 
is given by: 

/(s) = s + -. (41) 

s 

It is obvious that /(oo) = 00 and /(O) = 00 . Moreover, f'{s) = 1 — and 
for f'{s) = 0 we obtain the minimum value s = 1 which indicates an equal 
radii to the magnet and the current loop. For this case 6 = 2 + -^ and for 

the case that the current loop is put on the magnet 6 = 2 and g{2) = — §• 
For a ’’thick” magnet /ii —)■ 00 and thus also 61 —)■ 00 so that g{bi) —)■ 
which is unfortunately a rather slow decrease. Usually the magnet does not 
have to be too thick only enough to make \g{bi)\ < | 5 '( 6 )|. The function 
g{b) can be written explicitly in terms of the elliptic functions Ee{m) and 


^ _ h\ R2 Ri ^ ^ 

= irw + IT + IT- 

K1K2 -Til ii 2 


= \/ RiR 2 'i/bi — 2 cos 0'. 

(38) 



g(b) 



Figure 2: The function g{b). 


Ke{m) as: 

g(b) = l{i/632[-6£e(^) + (6 + 2)Jfe(^)] 

+ + {b - 2)Ke{-^-^)]} 

TT 

Ee{m) = / (16^/1 — msin^ 6 

Jo 

TT 

Keim) = f dO (42) 

Jo V 1 — m sin^ 6 

A graph of gih) is given in figure [2j 


5 Specific Examples 


Let us look at the case in which the radii of the magnet and the cnrrent 
loop are equal and the current loop is placed on top of a thick magnet. In 
this case: 


1 ^ 12,2 = 


IGtt R\ 


z. 


(43) 


Inserting equation (H3l) into equation (|2^ will result in: 



-iioMoh 


rRl. 


-Z. 


(44) 
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However, since the residual magnetic flux density is related to the magneti¬ 
zation as Mq = we may write: 


^ 9 .. R3 

Ft = --Brh^z. (45) 

o c 

We notice that for strong magnets of the NdFeB type the residual magnetic 
flux density ss 1 — 1.4 Tesla. Hence the pre factor is of order 1. 

Let us now look at the more general case of a thick magnet with a current 
loop on top but now the radii of the cylindrical magnet and the current loop 
are not equal hence h = 0, s ^ 1 and the force formula takes the form: 

^ 1 .. ??3 1 

Ft = --Brli^s^-^g{s -h -)z. (46) 

4 s 

Since: 

lim — -I —) = TTs (47) 

s—>oo S 

We arrive at the approximate force equation in the case that the radius of 
the magnet cylinder is much bigger than the radius of the loop current: 


TT 


Ft = -Brh 


. RIR 2 


z. 


(48) 


If we take the second derivative R to be of the order R 


we obtain: 


Ft ^ ^Brlp 



R 2 Z. 


(49) 


we see that the decisive factor is the ration of the current rising time and 
the time it will take a light signal to travel across the current loop. 

Finally we would like to address the question of the possibility of the 
device to lift from the ground for this the force generated by the device 
should be larger or equal to the gravitational force 

Fg = gm = gpV = gpLirRl (50) 

In the above g is the gravitational acceleration on earth (~ 9.8^ not to 
be confused with the function g{b)), m is the mass of the magnet, p is the 
mass density of the magnet and V is the volume of the magnet. The ratio 
of relativistic and gravitational forces is given by: 

R ^ -zBrR^s^-^g{s + ^) ^ BrR^ / -g{s + i) \ 

Fg gpLirs'^Rj AgpLn ^ ^ j ' 
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If the magnet is large with respect to the current loop we can take the limit: 


lim 

S —>oo 


[ ) 


TT 

72 


Which leads to the force ratio: 

Ft ^ 

Fg ~ AgpLs‘^ ■ 


(52) 


(53) 


Now since the current loops are small we may consider to put N loops on 

B? 

the magnet the maximum number in a single layer would be: N = -^ = 
hence the total force ratio is: 


NFt Brh^ 

Fg 4gpL 


(54) 


For a magnet of radius of 1 meter and thickens of 1 meter we obtain for a 
NdFeB magnet the density is p = 7500^ hence we need at least a ii ~ 
1.9 * 10^^^ for this Hying saucer to fly. This type of acceleration it typical 
to microwave currents of frequency 10 GHz and current amplitude of about 
4.8 Ampere. 


6 Conclusion 

We have shown in this paper that in general Newton’s third law is not 
compatible with the principles of special relativity and the total force on a 
system of a current loop and a permanent magnet system is not zero. Still 
momentum is conserved if one takes the field momentum into account. 

We have developed simple formulae for the total force in various cases 
including the case of equal radii of magnet and current loop and the case 
of a small current loop. This was extrapolated to the case of multiple loops 
and a specific example was done that the force of the system is equals its 
weight. 
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